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Abstract 
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I We construct weak axial one-boson exchange currents for the Bethe-Salpeter 

\^ I equation, starting from chiral Lagrangians of the NA(1236)7ryoaia; system. 

The currents fulfil the Ward-Takahashi identities and the matrix element of 

CN ■ 

^ ■ the full current between the two-body solutions of the Bethe-Salpeter equation 

CN , satisfies the PCAC constraint exactly. 
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I. INTRODUCTION 

The quantitative studies of the (anti)neutrino interaction with nuclei have now become 
reality, thanks to the existence of several neutrino detectors. Actually, the main goal of all 
^ , these studies is a firm empirical evidence for the neutrino oscillations |l|], 0, f^. 
^ ' Here we deal with the (anti)neutrino-deuteron reactions, 

u + d — > u + n + p , (1.1) 

: Ue + d — > e~ + p + p, (1.2) 

i7 + d — > V + n + p , (1.3) 

Ve + d — > + n + n. (1.4) 

In the recent experiment on the electron antineutrino-induced deuteron disintegration 
done at the Bugey reactor, the neutral current (NCA) reaction ( |1.3D and the charge current 
(CCA) deuteron disintegration ( |1.4| ) were probed at low energies (~ 1 MeV). While the 
measured cross section for the reaction NCA is in agreement with the previous measurements 
I^J^ and calculations 0J^, the new CCA cross section differs considerably from the previous 
results. 

Reactions ( |1.1| ) (NCN) and ( p..2| ) (CCN) are important for studying the solar electron 
neutrino oscillations. They are the main object of the SNO detector and should be 
measured soon. While the charge current effect will come presumably from the solar 
neutrinos {E,y < 15 MeV), the neutral current effect should reflect the interaction of all 
the neutrino species. If the ratio of the observed rates for reactions ( |1.2| ) and ( |1.1| ) would 
be measured significantly less than the standard model prediction, then the mixing of the 
neutrino species would be clearly demonstrated. 
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High energy neutrinos {E^, > 50 MeV) can be expected from solar bursts and from 
reactions of cosmic rays in the atmosphere (atmospheric neutrinos). 

The cross sections for reactions ( |l.l| )-( p~5D up to energies 170 MeV are presented in [jlO 



They were calculated within the frame of the standard nuclear physics calculations using the 
nuclear wave functions obtained by solving the Schroedinger equation with the commonly 
used nucleon-nucleon (NN) realistic potentials ||rT| . Besides the one-nucleon currents, static 
weak one-pion meson exchange currents (MECs) were also taken into account. 

At high energies, relativistic effects both in the potentials and currents (consistency of 
calculations) should be considered. In the standard non-relativistic approach, relativistic 
corrections are calculated in the 1/M expansion scheme (M is the nucleon mass). Actually, 
some of the components of the MECs (the time component of the vector MEC and the 
space component of the axial MEC) exist only as relativistic corrections of order 1/M^ 
which makes the calculations difficult. On the other hand, relativistic calculations can be 
carried out from the beginning within the framework of relativistic equations. This approach 
has recently been applied to the electron-deuteron interaction at high energies [Q, where 



the general discussion was done for the Bethe-Salpeter (BS) [|T^ equation and furthermore 
the electromagnetic interaction has been constructed in accord with the spectator (Gross) 
quasipotential equation. Here we apply the approach of the above mentioned paper to 
study the structure of the weak interaction in conjunction with the BS equation, having in 
mind to perform later a similar study for the Gross equation and to use the weak MECs in 
calculations of the neutrino-deuteron cross sections at high energies. 

The search for the MEC effects at low and intermediate energies has a long and successful 
history |I^-p!9|. The one-pion range MEC effects were firmly established in the reactions 
governed by the space part of the electromagnetic current and by the time component of 
the weak axial current which are of the same order ~ 1/M as the corresponding components 
of the one-nucleon current. Actually, the recent precise measurement of the transition 
rate for the reaction 

/i" + ^He — > u^ + ^H (1.5) 



and its comparison with the calculations |21] demonstrates clearly also the presence of the 



effect of the space component of the weak axial MEC. At higher energies, other meson 



exchanges p2| , p3| and relativistic effects start to play an important role. 

At low energies, the starting point for constructing the MEC operators are the low-energy 
theorems. In order to pass to higher energies, one employs chiral Lagrangians relevant to 
a nuclear system under investigation. In constructing the Lagrangians, symmetries play 
an essential role and dynamical symmetries imbedded in the Lagrangian reflect dynamical 
principles governing the system under investigation. Such a basic principle governing the 
physics of electroweak processes in nuclei at intermediate energies is the vector dominance 
model p3 . The vector meson flelds are introduced into Lagrangians flrst as massless Yang- 



Mills (YM) gauge flelds [^] belonging to the linear realizations of the chiral symmetry. 
Introduction of the mass terms for the vector meson flelds by hand makes this concept 
internally inconsistent. Also the interaction of the particles with the external flelds has to 
be introduced by hand since there are no other charges available which could be related to 
the external flelds. These defects of the concept of the massive YM flelds are later removed 
in the hidden local symmetries (HLS) approach [ P^ , p7| ]. 
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In this approach, a given global symmetry group Gg of a system Lagrangian is extended 
to a larger one by a local group Hi and the Higgs mechanism generates the mass terms for 
gauge fields of the local group in such a way that the local symmetry is conserved. For the 
chiral group Gg = [SU{2)l x SU{2)n]g and Hi = [SU{2)l x SU{2)r]i the gauge particles 
are identified [P§,P3] with the p- and ai mesons. An additional extension by the group U{l)i 



allows one to include the isoscalar u meson as well p8|. Moreover, external gauge fields, 
which are related to the electroweak interactions of the Standard Model, are included by 
gauging the global chiral symmetry group Gg. 

Let us note that the YM approach and the HLS concept with the vector meson fields 
belonging to the linear realization of the chiral symmetry are formally equivalent p6| , |27| . The 
linear and non-linear realizations of the HLS are related by the Stueckelberg transformation 
| 29| , |27[] . Lagrangians written in terms of fields of any realization are equivalent: the physical 
content of the model cannot depend on the parametrization of fields. 

Our starting point for constructing the weak axial MECs are chiral Lagrangians of the 
NA(1236)7rpaiC(j system. The first one is that with the YM vector meson fields belonging 



to the linear realization of the SU(2)xSU(2) chiral symmetry ||3T]]. It is a minimal La- 
grangian and contains no more than two derivatives of fields in each term. Another suitable 



Lagrangian constructed already within the concept of the HLS can be found in [^. These 
two Lagrangians differ by the vertex vrpai which brings into the MECs a model dependence. 
The model dependence of the electromagnetic isovector one-pion MEC upon this difference 
was studied in the backward deuteron electrodisintegration in It influences considerably 
the cross sections at energies about 800 MeV. 

The chosen Lagrangian has already been used in constructing the weak axial MEC oper- 
ator of the one-pion range within the extended S-matrix method |34]. In this method, 
the two-nucleon transition amplitude is represented by a set of tree Feynman graphs of the 
one-boson exchange type with the external nucleon lines on-shell. It is gauge- and Lorentz- 
invariant and it satisfies the PCAC equation. It is also assumed that the nuclear states are 
described by the Schroedinger equation solved with a one-boson exchange potential. The 
nuclear weak axial MEC operator is obtained after subtracting the first iterated Born ap- 
proximation from the two-nucleon transition amplitude. Let us stress that the first iterated 
Born approximation generally differs from the positive-frequency part of the nucleon Born 
term entering the two-nucleon transition amplitude in higher order in 1/M. As mentioned 
above, this MEC operator has been employed in |^ to calculate the weak axial MEC effects 
in reaction ( |1.5| ). 

Here we construct the weak axial MEC operator for the case when the two-nucleon 
nuclear states are described by the BS relativistic equation. In Sect.|T|, we present the 
Lagrangians. The transition operator is constructed in Sect.|Tl| starting from the YM La- 



grangian and its model dependence is studied in Sect.f^. For the construction, the same 



Feynman diagram technique is employed as in [^, but the nucleons are off-shell. In order 
to get the correct MEC operator, one should omit the nucleon Born terms from the transi- 
tion operator, which are already accounted for when calculating the matrix element of the 
one-nucleon current between the nuclear wave functions. We show that the full BS current 
satisfies the correct Ward-Takahashi identity (WTI) and that the resulting matrix element 
of the MEC operator between the two-body solutions of the BS equation satisfies the PCAC 
constraint exactly. 
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Similar program for the electromagnetic interaction has been done earlier in Ref. pS 
However, Gross and Riska use the minimal substitution in the kernel of the relativistic 
equation for generating the potential exchange currents, while the non-potential currents 
are derived from vertices. We derive both kinds of the currents from chiral Lagrangians. 
Generally, some difference can appear, since the currents derived from Lagrangians might 
be non-minimal. 

Our results and conclusions are given in Sect.|V]. 

The investigation made here is based on the normal Lagrangians. The anomalous La- 
grangians describe processes with the change of the internal parity of particles and they are 
related to the phenomenon of chiral anomaly I^U] . In our case, the vertices Tipuj and pu ai 
would play a role. They are characterized by the presence of the tensor e"^"'^ causing the 



37[] . However, 



change of the internal parity. The relevant Lagrangians can be found in 
the construction of these currents is beyond the scope of this paper. 

The weak axial MECs based on the Lagrangian approach have recently been studied in 
other models in p8| , p9 . 

We use the metrics defined by ( H ) . 



II. CHIRAL LAGRANGIANS 



We construct the MEG operators from a chiral Lagrangian of the Nirpaiuj system. Its 
essential part is the Lagrangian of the subsystem npai. The conventional Lagrangians of the 
Tcpai system contain trilinear and quadrilinear derivatives of the fields. They were studied 
already in sixties within the framework of the YM massless compensating fields. As it was 



shown in a minimal effective Lagrangian C^p^^ of the Trpoi system can be constructed so 



that it contains no more than two field derivatives in each term. This restriction determines 
the Lagrangian completely and it describes reasonably the related elementary processes up 
to 1 GeV. This construction was extended in ||3^, to the iVAvrpai system. 

However, the approach of the massless YM fields is internally inconsistent, as the masses 
of the heavy mesons should be put into the theory by hands and they violate the underlying 
symmetry. It was realized later |26|, that the internal inconsistency of the YM approach 
can be removed in the theory of HLS. In this theory, the masses of the heavy mesons are 
generated by the Higgs mechanism and they do not violate the symmetry any more. The 
HLS approach was applied also to the construction of the Lagrangian of the vrpai system 
| |26| | , where the high energy behaviour of the amplitudes was corrected by a different choice 
of correction Lagrangians in comparison with [|^, thus introducing a model dependence. 
However, the construction [^] was done for the heavy mesons on-shell and the Lagrangian 
cannot be used for constructing the MEG operators. Moreover, it also provides for the 
anomalous magnetic moment of the ai meson (5a„=+3, in contradiction to its commonly 
accepted value 5an = —1 H- Both gaps in construction have recently been removed . 



It has been also shown in , that the Lagrangians of the YM approach are formally 



equivalent to those constructed from the fields of the linear realizations of the HLS and that 
they are a sub-manifold of the manifold of the BKY Lagrangians. We first present the 
A^Avrpai Lagrangian C^^ of the YM type 0, 0, 



It can be easily extended to 
contain also the uj meson field. Then we introduce the HLS Lagrangian and 
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discuss possible source of the model dependence of our MEC operators. 

The YM Lagrangian of the NAnpai system can be written with an arbitrary mixing 



of the pseudoscalar and pseudovector vrA^A^ couplings |^ but we restrict ourselves to the 
pseudoscalar one only. With the u meson included, it is of the form 



fM _ r-M _|_ r-M , r-M 



npai ' 



(2.1) 



Nnp aiuj 



2/. 



9p,T, 



« 9 



pu 



TX 



f2.2) 



A/(0) 
NiTp aiuj 



^ (2 ^ - M) ^ - z^^TT ■ r75^ - c/AC/p^7^75(r ■ a^)^ 



Here we define 
P 



1 



pu Ppu + X + X dull -dyX df,n + n x 3,,^] + 0{'k'^) 

Jit 



and 

Ppu = d^py - dt^Pi,. 
Furthermore 



"pi/ 



p 1 



pv 



nM 

NAnpai 



^liEA^f^f^ ■ +igp^^''j,YT^ ■ Ppu + h.c. 



(2.3) 



(2.4) 



as 



Mt rrit 



Here the transition isospin operator T is defined 

(f) 

Further in Eq. ( p^ 





I . (2.5) 



(2.6) 



As it was discussed in |^0|, Gi ~ 2.6. Actually, another term ~ G2 present in the original 
Lagrangian was neglected in (|2.4|) because of one additional factor ~ 1/Ma. 



5 



The Lagrangian of the vrpoi system is of the form 

-p" X a'') ■ df,a^ - yP^" ' ^ ^yii + ^tt x a^^) . (2.7) 
In this approach, the hadron currents are mediated by mesons, 



ml „ -„ ml ^„ -.,„ ml 



jlr = — = — ^p^/^/ = — '-s^ - /.^^^ - x t? + o( 

9<^ 9p 9p 



TX 



(2.8) 

The Lagrangian ( |2.1D together with the currents ( |2.8| ) consist of all the vertices 
necessary for constructing the operators of the isovector vector and axial one-nucleon and 
TT, p, Oi and uo MECs in the tree approximation. 

The HLS Lagrangian consists of three terms like the Lagrangian C^^ , but now the 
heavy meson fields represent the nonlinear realization of the HLS for the groups Hi = 
[SU{2)l X SU{2)r]i X U{l)i. In comparison with (^ we have 

J^N.pa,^ = C-lfpa,. - 9p§j^ ^ {r ■ /V X n) v[/ + g^^A ^^.(^ ■ a, X U)^ 

(V \ ■ V 

I .!!£_ crM^JF(p) 1 . fiT/ + Iq A.^^^.o-Mi' ( -pip) . ft W 
^ 4M 4^^M2M^' ^ ^ 

+ 0{n'), (2.9) 

with 

^pJ = Ppy + ^P^M X ' (2-10) 

and the Lagrangian C^J^^^aiuj formally coincides with C^^fpa^u)- The same holds also for 
^^Anpai ■ As to the Lagrangian C^p^-^ , it is of the same form as given in Eq. (|2.7| ) with the 
last two terms changed to 

A£f,,^ = [p'. ■ Sp X d'^li + \pp ■ X a>t) . (2.11) 



The currents- analogues of Eqs. (|2.8| ) are given by the following set of equations 

771 ^ Tfi^ 

jSp ^ — ^^p^ — p_^p _ 2j^^^M X n + o(n2), (2.12) 

9u; 9p 

2 

j5p ^ p_^p _ j^gp^ _ ^j^g^^p ^ n 

9p 

+ - (-T^d'-U- gpS" + eA") x p'^, + 0{U^) . (2.13) 

In contrast to the YM currents Eq. ( p.8|) , the HLS curents contain also the external weak 
axial field A'^. 

Checking the Lagrangians and currents for the possible source of difference between the 
MECs one finds that it can arise from the vrpai and NNvrai vertices. 
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III. THE AXIAL MEC OPERATOR FOR THE BETHE-SALPETER EQUATION 



We first write the operator of tlie weak axial current for the ith nucleon (i= 1, 2) 



9a 



m 



l^A^,:^{q)P,:: + ^UAUq')q'tt■ 



Here is the weak axial coupling (the nucleon axial charge), g^ = -1.256, = p' i — Pi, 
the vector-meson propagator is generally designed as 



1 



and the pseudovector and pseudoscalar vertices are defined as 
The divergence of the current ( |3.1|) is 



(3.2) 



(3.3) 



^^m^A^(g2)f^ 



(3.4) 



Here the operator of the nucleon axial charge 6^(1) and the inverse of the nucleon propagator 
G~^{p) are defined as 



CAl^j = 9A I75y ) , G 



- M 



(3.5) 



The general structure of the weak axial MEC operators is given in Fig. 1. The weak axial 
interaction is mediated by the meson B which is either vr or ai meson. The range of the 
current is given by the meson B2 which is here vr, p, oi or u meson. In the graphs la and lb, 
TV is either for the nucleon N (in this case the graphs are the standard nucleon Born terms) 
or for the A(1236) isobar. We shall write the corresponding currents as Jg^(A/', B).ln the 
case of the BS equation, the nucleon Born terms do not enter the nuclear MECs because 
they are fully taken into account when calculating the matrix element of the one-nucleon 
current with the solutions of this equation. The graph Ic represents one type of contact 
terms which we shall denote as J"^^ [B) ■ It is connected with the weak production amplitude 
of the B2 meson on the nucleon. Another type of the contact terms is given by the graph 
Id where the weak axial current interacts directly with the mesons Bi and B2. We shall 
write it as JbiB2- The graph le is for a mesonic current which we shall write in the form 
^3132^^)- 'T^^ associated pion absorption amplitudes will be written as j(/l%^{N')^ ■^cB2 
and -Mbi_B2' ■^'^^ convenience, we shall modify the notations in some cases. 



A. The axial MEC operator of the pion range 



This operator is closely related to the one derived in 



from the Lagrangians C^^ 



Nnpaiuj 

and C^pai- Essentially, it has potential and non-potential parts, with the nucleon Born 



7 



terms in the potential part omitted. Let us note that the Lagrangian used for the derivation 
of the 71 MEC operator corresponds to the non-hnear chiral model with the pseudoscalar 
ttNN coupling which is connected to another model with any mixing of the pseudoscalar 
and pseudovector couplings by the unitary Foldy-Dyson transformation which makes all the 
models physically equivalent. 

In the given case, the potential amplitude consists of several contact terms where the 
weak interaction is mediated either by the ai- or vr meson. The ai-pole contact terms of the 
type J"^(ai) are 



with 



1 



4/. 



A;^;(g)e— r™(gOA^fe')r^ + (1^2) 



(3.6) 



(3.7) 



and 



A:^;(g)e— (7.r'")iA^(g2')f^ + (1 



2) 



The TT meson contact terms of the type J^^{it) are of the form 
where the pion absorption amplitudes are 



AM", 



and 



2 \mJ 
4 \mJ 



J^rM = iUq^Al{q^)^5anAUql)T'i + (1^2) = lU q^ Al{q^) M^^ 



/l2TTAUqi)Y^, 



2) 
2) 



(3.8) 



(3.9) 
(3.10) 



(3.11) 
(3.12) 

(3.13) 



In contrast to ref. [33|, besides the current (|3.9|) we have an additional term (|3.10 ) which 
disappears for the on-shell nucleons. 

A direct calculation of the divergence of the currents ( p.8|) , (|3.9| ), (|3.10|) and ( p.l3| ) yields 



9m 



7r 



2 A TT f ^2 



zUmlAUq'' 



M" +2AM'' 



f^q'Al{q')AM'::^^-zf^Ml^. 



It can be shown that 



-lUMl^ = e^(l) + e^(2), K 



+ 



(3.14) 



(3.15) 



where the one-pion exchange potential K- is defined as 

Taking into account Eqs. (|3.1CI|) and (|3.15|) , we can write Eq. ( p.l4|) in the form 



(3.16) 



C2 IT 



+ 



eA(l)+eA(2), K 



(3.17) 



This continuity equation allows us to define the potential axial tt MEC operator related to 
the BSE as 



(3.18) 



Let us note that the current (|3.6| ) does not enter the resulting potential current ( 3.18| ). As 
in the on-shell case it satisfies the continuity equation together with the non-potential 
currents (see below) which are constructed from the irpai vertices. So we find its natural 
place among these non-potential currents, which we will consider now. They are analogues 



of the amplitudes given in Sect. 4 of Ref. [Q. The explicit form of the currents is 



ml. 



8f: 

AUql)^ + (1 2), 
2^g^A^(g2) \m;^ + AM 



f u{q) ^'^™"r™ (gi) A^(gi2) I (gi ■ q, - ql)g^^ + ^^(g? - qV 



m 



m 



" e^^y,,T'(^{q,) 1 - 2(gi ■ g2)A^(g^) Al{qi)T:^ + (1^2 



P ^'^2^ 



rap 



m 



4/. 



^e'^"^^A^;{q,)TZ{qi)AUql)AUql)r^2 



2) 



(3.19) 
(3.20) 

(3.21) 
(3.22) 
(3.23) 



Strightforward calculation yields for the divergence of the currents ( |3.6| ) and ( |3.19| )- ( ^.23[ ) 

J,7^(aO + >^(aO + iT(7r) + + AJ.^^np)] = tUmlAUq' 



where an additional current AJp^{np) defined as 

1 



■' '■ p-K ' ■' '■ p-K 



Amnp) 



4/. 



g'^A^(g2)e"'"" /i^rl^ A^^qDV^, + (1 - 2) 



and a new pion absorption amplitude 

= ^ e'^"^- [l - (gi ■ q2)AUql)] AUql) f ^ + (1 - 2) 



(3.24) 



(3.25) 



(3.26) 



appear for the off-shell nucleons. 

It can be verified that the sum of the currents Jp^(7r) and AJ^^(np) is 



whereas the other contact terms satisfy 



M"" +M 

p-K ' pw 



(3.27) 



(3.28) 



as it should be, if the continuity equation ( p.24 ) should hold also for the off-shell nucleons. 
Evidently, the non-potential axial vr MEC operator is defined by Eq. ( p.24| ) as 



JbLM = J,7^(ai) + iT(ai) + >^(7r) + >^ + A>^(H . 



(3.29) 



As noted above, the natural place for the current ^citl*^!) among the non-potential ones. 

We now consider another contribution to the non-potential currents of the pion range 
coming from the A(1236) excitation currents. These currents appear because of vertices 
Oi A^A and vrA^A and they are of the form 



7, 



TTiVA 



jr(A, vr) = tUq^ AUq) MUA) , 

/ r \ 2 

/ JttNA 



Ml{A) =\il^\ q^or{A), 
\ J 

or (A) = [(T+)"TO:^'^(P) + {T+yT'S^/iQ)]^ g2AA^(g2)n +(1^2) 
Here the propagator of the A isobar is defined as 



1 



- Ma + is 



1 2 



^ 3 ' ' 3Ml 



2 p"/ + 



3Ma 



(p°7^ - ^7") 



(3.30) 
(3.31) 
(3.32) 
(3.33) 

(3.34) 



and the transition isospin operator T is defined in accord with Eq. (|2.5|) . The total A isobar 
current of the pion range is 



with the divergence 

g,jr(A) = iUmlAUq)MUA) 
We define the BS weak axial vr MEC operator as 



J 



a p 
BStt 



ex 



JTsM + J%'sAnp) + ir(A) 



(3.35) 



(3.36) 



(3.37) 



with JbsAp)^ JbsMp) and J^''(A) given in Eqs. (|3l8|) , (ICTI) and (|05|) , respectively. 
Using Eqs.(|5T7|), (^) and (|CT| ) we obtain the WTI for the tt MEC operator 
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where the pion absorption amphtude A^"(2) is 



2AM'' + + +M 



(3.38) 



(3.39) 



At low energies, the weak axial potential vr MECs play little role. The non-potential 
current ( p. 231 ) dominates the time component of the weak axial vr MEC and its effect can be 
clearly seen in light nuclei Jl^, [|^- [|1^]. On the other side, the current (|3.30|) contributes 



significantly to the space component of this current, as observed in the weak reactions in 



lightest nuclei, particularly in the reaction ( |1.5| ) p0|j21 



At higher energies, as in the case of the electromagnetic interaction, heavier mesons are 
expected to play a non-negligible role. Now we consider the weak axial MECs of the p, ai 
and u range. 



B. The axial MEC operator of the p and ai meson range 

As it will become clear soon the p and ai exchanges should be considered in the given 
model simultaneously. The weak axial MECs of the p meson range derived from our model 
Lagrangian C follows. The only contact term is 



(3.40) 



where 



9 a 9p 



U V2 



^5„„(75<)ig2.A^^(g2)r^,(g2) 



(3.41) 



The mesonic currents are 



aTp(«i) = ^-7r9pm'p 



:9a ^2^2 
2 



(g. + gi.)A„f (g) A,^^(gi) - q^A,l'Jq)A^^{q,) 
-gicA.f (9)V,.(gi)] e^^-tl^Ap{q2)t-^{q2) + {1 ^ 2) 



(3.42) 



and 



(3.43) 



with 



M:^p{^) = - 



9a (9p 



(2g2C + gic)Aa^.(?i) - (2g2. + gi.)A,^^^(gi) 



e-^-tl^q^<Ap{q2)t"2,iq2) + {1 ^2). 



(3.44) 



Here we modify the notation for the pion absorption amplitude from A^^ip ■^aipi'^)- 
The operator for the p meson exchange potential is 
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9p 



rL(-g2)A;'^fe)fL(g2) 



For the current of Eq. ( p.40|) , we can write the divergence in the form 



(3.45) 



(3.46) 



Using the exphcit form of the amphtude A^^p from Eq. ( p. 411) , the second term in Eq. ( p. 46]) 
can be cast into the form 



an 
P 



(3.47) 



where V^"" is the part of the potential Vp given by the anomalous part of the pNN vertex 
entering the anticommutator of the potential with the nucleon axial charge. We denote the 
missing part of the anticommutator by — Ap. It is of the form 



6^(1) + e^(2), , =igA(^-^] e'^^^v\^Ar{q2)tiM + (1 - 2). 



A. 



In the next step of deriving the WTI for the p meson MECs we calculate 



(3.48) 



^Ap{q2)f%{q^) + (1^2)= if.Ml^pi 



X] . 



(3.49) 



However, the amplitude AA'^^p{x) does not coincide with the amplitude Maip{'^) from 
Eq. ( p.44| ) and we get 



0.1 p 



TT 



tUmlAl{q')Ml^p{7,) + [e^(l) + e^(2) , 



(3.50) 



"^o, p, 



which together with Eqs. (|3.46| ) and ( p.47| ) provides for the exchange current J bs p{^^)i 
defined as 



' BSp 

with the WTI in the form 



(3.51) 



1.0, p 
IpJ BS pi^^) 



e^(l) + e^(2), \/J + lUml AUq')Mp{2) + if^ (Ml,p{x) - M^^ ^ 



TT, 



(3.52) 



where the pion absorption amplitude is given as 



Mp{2) = m:p + Ml,p{7^). 



(3.53) 
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Using Eqs. (|3.44j) and ( |3.49|) , we derive for the difference of the amphtudes Maipi^) 
7W^^p(7r) the following equation 



TT 



(3.54) 



Evidently, this difference is an operator of the ai meson range. As we shall see below, it will 
be compensated by the contribution from the weak axial exchange currents of the ai meson 
range. 

Considered together with the p meson MECs, the oi meson MECs contain only 3 contact 
terms J,7,r7r) (i=l,2,3) 



The associated pion absorption amplitudes are 



9a 



In 

2 

^|^'^"^"(7M^")lAaTfe)f|; +(1^2), 



(3.55) 

(3.56) 
(3.57) 
(3.58) 



Let us note that the currents ( 3.55| ) are derived from the same NNvrai vertices as the tt 
MECs ( p.6| )-( |3l8D , only the role of the weakly and strongly interacting mesons is exchanged. 
The Oi meson exchange potential is 



5 n 



(3.59) 



Let us now calculate the divergence of the currents J".^^ (vr). In the first step we compute 

%j:raM) = iUmi/v^{q')Mi^,^ - ium:^^^ . (3.60) 

Using the explicit form ( p.56| ) of -Mciai follows that 



-^A-^:.a. = [eA(l) +e^(2),V;, 

We next calculate 



(3.61) 



C3 ai 



(3.62) 



Employing Eqs. ( |3.60|) - (|3.62|) , we have the following WTI for the weak axial Oi MEC derived 
so far 



qj7saM^) = [eA(l) + 6^(2) , Va,] +tUmlAUq')Ml{2)-tU {Mt^a^+M% 



(3.63) 
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^1 fl 



with the current JBSaii^^) defined as 



ai 



(3.64) 



i=l 



and with the pion absorption amphtude Ai (2) given by 



j=l 

The sum of the last two terms at the right hand side of Eq. ( p. 631) yields 



(3.65) 



ff;::A,7(gi)fL(g2) + (1 2) 



(3.66) 



Observing that the right hand sides of Eqs. (|3.54| ) and (|3.66|) are of the same form but with 
the opposite sign we derive the WTI for the sum of the weak axial p and ai MECs 



BSai 



6^(1) + 6^(2) , V, + Ki 



Using Eqs. ( |3.54| ) and ( p. 66] ) once more, we get from Eqs. ( |3.53| ) and ( |3.65D 



^,(2) + ^„^(2) = m% + m:,^ 



ci ai 



(3.67) 



(3.68) 



We now derive the A excitation currents of the p and ai range. In analogy with 
Eqs.( |3.3(]| )-( |3.36| ) we have for the p MEG the following set of equations 



m,r Ml ' 
{T^rT'^j,j^Sly{P) + (T+)»T"^);^(Q)757^]^ (q,,g 

-Q^9x') ^Uq2)TUq2) + (1 - 2). 



r/iA, ai) 
r/{A, ix) 
M%A) 
67(A) 



The total A isobar current of the p meson range is 

J7(A) = J7(A, ai) + J7(A, tt) , 

with the divergence 

qJ;^{A) = iUmlAl{q^)M%A). 
Similarly we have for the A excitation currents of the ai meson range 



(3.69) 
(3.70) 
(3.71) 

(3.72) 
(3.73) 
(3.74) 
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NA 



J,7(A, tt) = iUq^^l{q)Ml^{^) 



NA 



0,7(A) = [(T+rT'^S'/iP) + {T+TT^S'^/iQ)] Al,iq,)Tl; + (1^2). 
The total A isobar current of the ai meson range is 

rj{A) = rj{A, a,) + J:I\A, vr) , 

with the divergence 

g.^a?(A) = ^/.m^A^(g^)^^^(A). 
The total BS weak axial p+ai MEC is defined as 



JbS p+ai i^^) 



and it satisfies the WTI 

with the pion absorption amplitude ■M'^+aii'^) given by 

M^,^,,{2) = M% + Ml^ix) + M:^,^ + M;{A) + MliA). 
The last MEC operator we will now construct is that of the u meson range. 



(3.75) 
(3.76) 
(3.77) 
(3.78) 

(3.79) 

(3.80) 

(3.81) 
(3.82) 

(3.83) 



C. The axial MEC operator of the to meson range 

The u meson is a vector meson with the isospin zero and the only MEC operator which 
appears is the contact current JcHij^) 



Mt., = ^ ( ^) (75<g2.r")^ /\!^\q,)f,,{q,) + (1^2) 



2M 



V2 



(3.84) 
(3.85) 



Here f 2,^(52) is given by Eq. ( ^.7| ) with the isospin operator omitted and with the change 



The uj exchange potential reads 



9u 



V. = (^f j ri,(-g2)Ar(g2)r2.(g2). 

Calculating the divergence of the current ( ^.84| ) yields 



(3.86) 
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Q,J:!^{^) = iUm'7^^l{q')M:^ - iUM:^. (3.87) 



Using Eq. (|3.85|) for the amplitude A^"^ and Eq. ( |3.86|) for the oj exchange potential K 



we verify that the second term on the right hand side of Eq. ( |3.87| ) provides the needed 
anticommutator of the nucleon axial charge and the potential. Finally we have the WTI for 
the BS weak axial uj MEC 

= tUm'nAUq')M:^ + [6^(1) + 6^(2), K.]^ . (3.88) 

We have completed the derivation of the BS weak axial one-body operator and of the it, 
p, ai and u MEC operators in the tree approximation. These operators satisfy separately 
the WTI (p:4|),(|338|), and ( pISHQ . Starting from them, we verify in the next step that 

the matrix element of the total current satisfies the standard PCAC equation . 

D. The continuity equation for the current matrix element 

We define, in accord with Sect. 2 of Ref. the full BS weak axial current as 

Ji's = ^J^aA)G2' + <(l,2)Gr^ + Js'siex) = + J^^(ex) , (3.89) 

where 

r^Uex) = r^Liex) + J7sp+a,{ex) + J,l^(7r) . (3.90) 

Using the WTIs for the one- and two-nucleon currents yields for the divergence of the full 
BS current 

%J7s = [eA(l) + e^(2), g-^]+ + iUmlAl{q^)M\ (3.91) 
where the inverse Green function is 

g~^ = G-B\ + V, (3.92) 

with 

V = V. + Vp + Va,+V.. (3.93) 

The BS propagator in term of the single-particle propagators reads Gbs = —iGiG2 and 

= tTlG^' + t^G^'+M:i2)+M;^,^i2)+M:^. (3.94) 

Because the two-body BS wave functions for both bound and scattering states satisfy the 
equation 

g-^\^lj> = <iP\g-^ = 0, (3.95) 

the matrix element of the divergence of the BS current ( p.91| ) fulfil the standard PCAC 
constraint 

q, < V^l J^g \ij>= lU mlAUq') < > . (3.96) 
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E. The strong form factors 

We have deah with the point BNN vertices so far. The usual way to introduce the strong 
interaction effects in the BNN vertices is to introduce the form factors Fbnn{<^) with the 
normahzation FB7VAf(iiin)=l. Then the vertex defined in (p.3|) becomes 



BJ 

T^F^Mq')- (3.97) 
It is a simple matter to verify that the one-nucleon current ( p.l|) should be redefined similarly 

J-^(l^^)F^MQi)■ (3.98) 



Then for the divergence of this current we have the same Eq. ( p.4[ ) but with the vertex 
from ( |3.97|) and with the nucleon axial charge eA{i) redefined according to the prescription 

6^(0 — ^ e^ii) F^r^^iq^) . (3.99) 

We could also leave the axial charge unchanged but then we should redefine the propagator 
in a similar manner 

G-\p) ^ G-\p) F^NNiq-) =^-M-J:, (3.100) 

with 

S = (i.-M)[l-F^^^(gf)]. (3.101) 

Introducing the strong form factors into the potential currents ( |378|) and ( p.9|) -( prT3D and 
to the A excitation currents ( |3.30| )-( p.31| ) follows the standard prescription [^] and the pion 
propagator connecting the baryon lines is modified as 

AUq^) AUq^) Fl^^ii) . (3.102) 

This prescription is not suitable for the non-potential currents (|3.19|) - (|3.23|) . In this case the 
correct procedure is 

A|(g2) Af (g2) FsNNiq^) , 5 = vr, p . (3.103) 

In order that the WTI ( p.24| ) would be satisfied in the presence of the strong form factors, 
the pion propagator of the current ( ^^ should be modified as 



A-(g2) ^ A-(gf)F^;vjv(gD^piViv(g|) . (3.104) 

Checking the WTI for operators of the p and ai meson range we conclude that it would be 
valid also with the strong form factors included provided 



1. In all operators of the p range 



K{q}) — a^f(9- ) Fl^N{q-) ■ (3.105) 
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2. In the amplitudes A^cjai ^"^^ -^cgai 

pNN 



) — A- (q^) F'{q^) . (3.106) 



In the amplitude Mciai the potential Vai and in the currents J^l^i^A, ai) and 

^ A- (q^) F'M) ■ (3.107) 



Finally it holds for the operators of the ui meson range that the necessary modification is 

A^k') ^ ^FiQ-) F^NNiQ-) ■ (3.108) 

In the next section, we investigate the model dependence of our MECs appearing due to 
the difference in the physical content of the YM type minimal Lagrangian C^^ and of the 
HLS Lagrangian £^ . 



IV. MODEL DEPENDENCE OF THE BS MECS 



The derivation of the MECs from the HLS Lagrangians and currents ( |2.9| )- ( |2.13 ) is 



analogous to that of the previous sections. Therefore, we present here briefly the results and 
emphasize the difference leading to the model dependence of the currents. The currents of 
this section will be labeled by an additional label H in order to distinguish them from the 
currents derived from the minimal Lagrangian 



A. The axial MEC operator of the pion range 

It is strightforward to check that the potential axial n MEC operator JBSniP^ coincides 
with the current Jbs7t{p) ol Eq. (|3.18| ). But generally, the non-potential currents differ. It 
holds only for the new pnir current Jp:j^{iT, H) that 

r,^{7T,H) = j;^{7t), (4.1) 

where J^^(7r) is given in Eq. ( |3.2CI| ). Checking further the currents J^/ of Eq. ( |2.8| ) and J^'^ 
of Eq. (|2.13|) we find that the new pvr contact term is now 

j;i^{l,H) = 2J;^^ (4.2) 

where the current is given in Eq. ( p.23|) . Here we label the current (|4.2|) additionally, 
because another pn term appears due to the presence of the vertex ~ d'^U x p ^ in the current 

j;,^(2, H) =^ e'^^^TTM K 92 - (Qi ■ 12)9^ A^(g2')f ^ + (1 2) . (4.3) 
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The new current j°^(ai, H), obtained using the Lagrangian A£^^^ of Eq. is 



+ ^ A^Hg2) (gi)A;,(gO [q'^q^ - (q ■ q,)g^^] ^''M)^2 + (1 2) , (4.4) 



The new non-potential currents satisfy the following WTIs 







(4.5) 



J,T(1, H) + j;^^(vr, H) + ^J;^{np) = mlA^q 



M" +M 



(4.6) 



where the current AJ^^(np) is given in Eq. ( p.25| ) and the pion absorption amplitudes A4p^ 
and Aip^ are from Eqs. ( |3.21| ) and ( |3.26| ), respectively. In analogy with the non-potential 



currents of the Sect. [Ill A , in addition to Eq. ( p.27|) we have now 

q,j;i:il,H) = lU 

instead of Eq. ( ^ ). 



A^" +M 

■' '■p-K ^ ■' '■pn 



(4.7) 



B. The axial MEC operator of the p and ai meson range 

In this model, the p and ai axial MECs satisfy the WTI separately. We first consider 
the p meson axial MECs. 

Instead of the current Jcpij^) of Eq. ( p.4(]| ), we have now two contact terms J^^piji^H) 
and4';(7r,i7) 



cp 



and 



where the pion absorption amplitude Xi'^^p{H) is 

M^c.piH) = -^j-[^ys^--Tl^A;^iq2)TUq2) + (1 - 2) 



(4.8) 



(4.9) 



(4.10) 



The first contact term appears due to the piece ~ p'^ x 11 in the current ( |2.13D , while the 
second one comes from the second term on the right hand side of the Lagrangian ( p. 91) . It 
also holds that 



(4.11) 
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with the current ja^p{a,i) given in Eq. ( p.42| ). An additional contact current is now generated 
from the vertex ^ a" x in the current ( p.l3|) 



J^l^^iH) = ^gAe'^"^-A:^{q,)tl^ U A^^^iq,) - q.^A^^^iq,)] f^^fe) + (1 2) . (4.12) 



The last current we consider is the mesonic current 

JZi^,H) = zUq^Al{q')M:^^{H) 

with 

9a (gp 



+ {-2q'2 - gi ■ 92 + qt) K,\qi) \ ^;\q2) r^,(g2) + (1 2 

It is strightforward to find that the the p meson axial MECs of this section 

^1 M 



satisfies the WTI 



J55^(ex, H) = [eA{l) + e^(2), Vp\ ^ + lUmi A^(g^)A^^(2, H) 
with the pion absorption amplitude given by 

m]{2,H) = M:^p{H) + M:^p{H) + Ml^piH). 
As to the Oi meson axial MECs, we have now only the current 

jZ,^iex,H) = JZ.in), 
of the currents (|3.55| ), for which the WTI ( p.60| ) is valid. 



(4.13) 



(4.14) 



J^spiex, H) = J,7^(vr, H) + J.^/vr, H) + J,7^(ai, H) + J^iH) + J^^^^in, H) , (4.15) 



(4.16) 



(4.17) 



(4.18) 



C. The continuity equation for the current matrix element 

In analogy with Eq. ( p.89| ), the full BS weak axial current is now 



where 



Js^{ex,H) 



JZAex,H) + rB'spiex,H) + J^sa^ex^H) + J^i 



Jssa^ex^H) + J,7(A). 



(4.19) 

(4.20) 
(4.21) 
(4.22) 
(4.23) 
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The associated pion absorption amplitudes read 

M\H) = zf-^Ga 1 + zf^Cr' + Ml{2) + A^^(2, H) + M^p, H) + M^^ , (4.24) 

M%2, H) = MI{2, H) + M%/\) , (4.25) 
Ml^{2,H) = M:^^^ +Ml{A). (4.26) 

An analogue of Eq. (p.91|) is then 

q.rB^siH) = [e^(l) + e^(2), Q-'^ + tUml^Hq^) M\H) , (4.27) 

which again leads to the standard PCAC constraint 

q,<i^\r/s{H)\i^>= lUml/VM^) < >l'^(/f) |^ > . (4.28) 

The possible model dependence is given by effects due to the difference between the 
currents J%g{ex) and J'^^^ex, H), which reduces to investigating the difference between the 
currents 

r^{M) = JZiai) + J;i:{ai) + JZi^) + JZ. (vr) + J^^ (tt) , (4.29) 

and 

r^iH) = 1 j;^^(l, H) + r^^{2, H) + j;^^(ai, H) + J^^iH) + J,7^(7r, H) + J,7^(7r, H) . 

(4.30) 

The presence of the vr meson axial MECs entering these two equations may indicate that the 
model dependence can appear already at relatively low energies. However, inspecting the 
structure of the currents Jci'^r('^i) l^pni^y-^) shows that the leading terms cancel each 
other and only the momentum dependent terms survive. So the model dependence should 
be expected presumably of the short range nature and can be expected at high energies. 



V. RESULTS AND CONCLUSIONS 

In this paper, we have investigated the structure of the axial MECs for the NAirpaiUJ 
system in the conjunction with the BS equation. Our starting points are Lagrangians and 
one-body currents of this system (Sect.|I|). One of the Lagrangians, is the minimal one 
and it is approximately invariant under the local [SU{2)l x 5'f/(2)/j] x f/(l) symmetry. The 
heavy meson fields are introduced as massless Yang-Mills compensating fields belonging to 
the linear realization of the chiral symmetry. The symmetry is violated by the heavy meson 
masses, which are introduced, together with the external electroweak interaction, by hands. 

Another choice considered here, is the Lagrangian , which reflects the [SU(2)l x 
SU{2)p']i X U{l)i hidden local symmetry. In this case, the symmetry is not violated by the 
heavy meson mass terms and the external electroweak interaction enters naturally into the 
scheme as gauge fields of the initial global chiral group [SU{2)l x SU{2)ji]g. We take the 
heavy meson fields belonging to the non- linear realization of the HLS, having in mind that 
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the Lagrangian expressed in terms of the fields belonging to the linear realization of the 
HLS is physically equivalent. Let us remind that the transition between the representations 
is done using the Stueckelberg transformation. 

Essentially, both Lagrangians differ by the choice of correction Lagrangians needed to 
obtain correct amplitudes of elementary processes up to the energies ~ 1 GeV. 

In Sect.|T|, (a) the tt, p, ai and uj meson axial MECs from the minimal Lagrangian 
are derived. These currents satisfy the WTIs separately but the case of the p and ai meson 
axial MECs, which should be considered together, (b) the WTI for the full BS current and 
the divergence of its matrix element between the two-body BS wave functions is given and 
it is shown that it satisfies the standard PCAC constraint and (c) the introduction of the 
strong form factors into the BNN vertices is discussed. 

The derivation of the analogous currents for the Lagrangian is presented in Sect. [TVI, 
where also the model dependence is indicated. 

The derived currents, with the nucleon Born terms added and sandwiched between the 
Dirac spinors, can be used in the standard nuclear physics calculations. To our knowledge, 
more complete set of the axial MECs has not yet been published. 



In comparison with earlier works , , we can say that the set of the axial MECs 
derived here is more realistic. Indeed, besides the ai meson exchange, the considered tt, p 
and UJ meson exchanges enter the realistic NN potentials of the OBE type [53|. Let us note 
that the ai meson exchange was included into OBEPQB in [0]. From the meson exchanges 
which are still not included, only the a meson exchange is of importance for OBEPs in 
describing the NN attraction at medium distances. This particle is controversial and it is 
difficult to include it consistently. It enters the linear a model naturally as a partner of 
the TT meson. However, the mass of the resonance |45] having the quantum numbers of the 
cr meson is by a factor ~ 2 larger than its mass extracted from the fit of OBEPs to the 
data On the other hand, in the non-linear a models, to which our models belong, this 
particle does not enter by construction and can be put in by hands only. 



Besides, Ref. contains an error |^ in sections dealing with the non-linear a meson 
model. Indeed, it should be {p'^ — Pi — 2g)y in Eq. (40) instead of {p'l — pi — q)u, which 
removes (a) the confusion about the axial current nonconservation even in the chiral limit 
given in the paragraph after this equation and (b) the q dependence on the right hand side 
of Eq. (60). It is seen then that the correct current can be simply redefined so that it satisfies 
the WTI and the matrix element of its divergence PCAC. 

As the next step would be calculations of the cross sections for the weak reactions like 
Eqs. (prT|)- ([L^ ). In order to do them consistently, one has to use solutions of the BS equation 
for the nuclear wave functions using a OBEP constructed in the same model. Such calcu- 
lations would give more confidence in conclusions about the charged and neutral processes 
in neutrino disintegration of the deuteron and consequently, about the neutrino oscillations. 
These are important considerations for understanding the solar neutrino experiments at 
SNO and other possible experiments at higher energies. 
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Fig. 1. The general structure of the weak axial MEC operators considered in this paper. The 
weak axial interaction is mediated by the meson B which is either tt or ai meson. The range 
of the current is given by the meson B2 which is here tt, p, Oi or uj meson. The graphs a, b, 
represent the current jy^{N' ^ B) with M either for the nucleon N or for the A(1236) isobar. 

The graph c represents a contact current J^^^{B). It is connected with the weak production 
amplitude of the B2 meson on the nucleon. Another type of the contact terms is given by the 
graph d, J^fsj, where the weak axial current interacts directly with the mesons Bi and B2. 
The graph e is for a mesonic current jy^Q^{B). The associated pion absorption amplitudes 
correspond to the graphs where the weak axial interaction is mediated by the pion, but with 
the weak interaction wavy line removed. There are three types of these amplitudes in our 
models: M%^{M), M^b^ and M%^b2- 
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